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ABSTRACT 



' A comparison is presented of Sunyaev-Zeldovich measurements for 1 1 galaxy clusters as obtained by Planck and by the ground-based interferom- 
eter, the Arcminute Microkelvin Imager. Assuming a universal spherically-symmetric Generalised Navarro, Rrenk & White (GNRW) model for the 
" cluster gas pressure profile, we jointly constrain the integrated Compton-K parameter (Y^oo) and the scale radius (6500) of each cluster. Our resulting 
, constraints in the K500 - ^soo 2D parameter space derived from the two instruments overlap significantly for eight of the clusters, although, overall, 
fS| ■ there is a tendency for AMI to find the Sunyaev-Zeldovich signal to be smaller in angular size and fainter than Planck. Significant discrepancies 
exist for the three remaining clusters in the sample, namely A1413, A1914, and the newly-discovered Planck cluster PLCKRSZ G 139.59-1-24. 18. 
The robustness of the analysis of both the Planck and AMI data is demonstrated through the use of detailed simulations, which also discount 
confusion from residual point (radio) sources and from diffuse astrophysical foregrounds as possible explanations for the discrepancies found. 
>v>( ' Ror a subset of our cluster sample, we have investigated the dependence of our results on the assumed pressure profile by repeating the analysis 
' adopting the best-fitting GNRW profile shape which best matches X-ray observations. Adopting the best-fitting profile shape from the X-ray data 
\ does not, in general, resolve the discrepancies found in this subset of five clusters. Though based on a small sample, our results suggest that the 
- ■ ■ adopted GNRW model may not be sufficiently flexible to describe clusters universally. 

Key words. Cosmology: observations - Galaxies: cluster: general - Galaxies: clusters: intracluster medium - Cosmic background radiation. 
X-rays: galaxies: clusters 



1. Introduction 

Clusters of galaxies are the most massive gravitationally bound 
objects in the Universe and as such are critical tracers of the for- 
mation of large-scale structure. The size and formation history 
of massive clusters is such that the ratio of cluster gas mass to 
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total mass is expected to be representative of the universal ra- 
tio, once the relatively small amount of baryonic matter in the 
cluster galaxies is taken into account (e.g.. White et al. 1993). 
Moreover, the comoving number density of clusters as a func- 
tion of mass and redshift is expected to be particularly sensitive 
to the cosmological parameters cr^ and Q.^ (e.g., Battye & Weller 
2003). 
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The Sunyaev-Zeldovich (SZ) effect (see Birkinshaw 1999; 
Carlstrom et al. 2002 for reviews) produces secondary 
anisotropics in the cosmic microwave background (CMB) ra- 
diation through inverse-Compton scattering from the electrons 
in the hot intracluster gas (which also radiates via thermal 
Bremsstrahlung in the X-ray waveband) and the transfer of 
some of the energy of the electrons to the low-energy photons. 
Moreover, the surface brightness of an SZ signal does not depend 
on the redshift z of the cluster Hence an SZ-effect flux-density- 
limited survey can provide a complete catalogue of galaxy clus- 
ters above a limiting mass (e.g., Bartlett & Silk 1994, Kneissl 
et al. 2001, Kosowsky 2003, Ruhl et al. 2004). 

Analyses of observations of galaxy clusters via their SZ ef- 
fect. X-ray emission or gravitational lensing are often based 
on some spherically-symmetric cluster model in which one as- 
sumes parameterised functional forms for the radial distribution 
of some cluster properties, such as electron density and tem- 
perature (Sanderson et al. 2003; Vikhlinin et al. 2005, 2006; 
LaRoque et al. 2006; Feroz et al. 2009b; AMI Consortium: 
Zwart et al. 2011; AMI Consortium: Rodriguez-Gonzalvez 
et al. 2011; AMI Consortium: Hurley-Walker et al. 2011; AMI 
Consortium: Shimwell et al. 2010), electron pressure and den- 
sity (Nagai et al. 2007; Mroczkowski et al. 2009; Arnaud et al. 
2010; Plagge et al. 2010; Planck Coflaboration VIII 2011), or 
electron pressure and entropy (AMI Consortium: Olamaie et al. 
2010; Allison et al. 2011). 

The motivation for this paper is to augment SZ measure- 
ments obtained with Planck ' for a sample of 1 1 galaxy clus- 
ters with refined higher-resolution SZ measurements obtained 
with the Arcminute Microkelvin Imager (AMI) interferometer 
Such a combination is an interesting and potentially very power- 
ful way to pin down the gas pressure profile of individual galaxy 
clusters as it relies on a single well-understood astrophysical ef- 
fect. In addition, Planck and AMI SZ measurements exploit very 
different aspects of the SZ signature: Planck eff'ectively uses its 
wide frequency coverage to identify the characteristic frequency 
spectrum of the SZ effect while AMI exploits its higher angular 
resolution to perform spatial filtering to identify SZ clusters and 
constrain their parameters. Combining measurements by these 
two instruments not only provides a powerful consistency check 
on both sets of observations, but may also break, or at least re- 
duce, the observed parameter degeneracy between the derived 
SZ Compton-y parameter and the cluster angular size which of- 
ten results due to the finite resolution of SZ telescopes (Planck 
Collaboration VIII 2011). 

The paper is organised as follows. In Section 2, we out- 
line how we selected our sample of 1 1 galaxy clusters for this 
comparison work. In Sections 3 and 4, we describe the Planck 
and AMI observations of our cluster sample, respectively. In 
Section 5 we present the pressure profile that we have used to 
model the clusters and constrain parameters. The analysis of the 
real data from both experiments is also described in Section 5. 
We follow this by validating our analysis methodology and in- 
vestigating the effects of diffuse foreground emission on the 
Planck constraints with simulations in Section 6.1. Section 6.2 
presents a similar simulations-based investigation of the effects 
of residual point sources and analysis methodology on the con- 



' Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions form NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 



straints derived from the AMI interferometric data. With a view 
to explaining some of the discrepancies we observe, in Section 7, 
we investigate the possibility of relaxing the assumptions regard- 
ing the universal pressure profile adopted for our cluster sample. 
Here we also examine the consistency of the Planck and AMI 
SZ results with complementary constraints from high-quality X- 
ray observations for a subset of our cluster sample. We conclude 
with a discussion in Section 8. 

2. Selection of the cluster sample 

An original sample of 26 clusters was defined at the beginning 
of this study. 24 of the clusters were identified as members of 
the sample by virtue of the fact that they were both present in 
the Planck Early Sunyaev-Zeldovich (ESZ) cluster catalogue 
(Planck Collaboration VIII 201 1), and had also already been ob- 
served and detected with AMI during the course of its normal 
observing programme. Note that these 24 clusters had been ob- 
served by AMI as part of differing scientific programmes and 
while each programme had a well-defined sample, the result- 
ing set of clusters used in this paper does not constitute a well- 
defined or complete sample. 

To this sample of 24, two newly-discovered Planck clus- 
ters were added, for which AMI made follow-up observations. 
The complete list of the original cluster sample, their coordi- 
nates and redshifts is presented in Table 1 . The sample was then 
screened to include only clusters that had (i) a firm SZ detec- 
tion by AMI (> 3cr before source subtraction and > 5cr af- 
ter source subtraction) and (ii) a benign environment in terms 
of radio point sources (we discard clusters with total inte- 
grated source flux densities greater than 5 mJy within a radius 
of 3' or greater than lOmJy within a radius of 10' from the 
phase centre). This reduced the sample to 11 clusters span- 
ning a wide range in redshift, 0.11 < z < 0.55 . The clusters 
in the new sample are A2034, A1413, A990, A2409, A1914, 
A2218, A773, MACS Jl 149+2223, RXJ0748-I-5941, PLCKESZ 
G139.59+24.I8, and PLCKESZ G121. 1 1+57.01. The sample 
includes two cool-core clusters, A1413 and A2034 (Pratt & 
Arnaud 2002; Kempner et al. 2003; Vikhlinin et al. 2005; 
Govoni et al. 2009) and two newly-discovered Planck clus- 
ters, PLCKESZ G139.59+24. 18 and PLCKESZ G121.1 1+57.01 
(Planck Collaboration VIII 2011). The last two have been ob- 
servered in the optical with the RTT-150 telescope- as part of 
the Planck follow-up programme. The resulting spectroscopic 
redshifts measured for the brightest cluster galaxies within these 
two clusters are given in Table 1 . 

3. Description of Planck data 

Planck (Tauber et al. 2010; Planck Collaboration I 201 1) is the 
third generation space mission to measure the anisotropy of the 
cosmic microwave background (CMB). It observes the sky in 
nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular resolution from 31' to 5'. The Low Frequency 
Instrument (LFI; Mandolesi et al. 2010; BersanelH et al. 2010; 
Mennefla et al. 201 1) covers the 30, 44, and 70 GHz bands with 
amplifiers cooled to 20 K. The High Frequency Instrument (HFI; 
Lamarre et al. 2010; Planck HFI Core Team 2011a) covers the 
100, 143, 217, 353, 545, and 857 GHz bands with bolometers 
cooled to 0. 1 K. Polarisation is measured in all but the highest 
two bands (Leahy et al. 2010; Rosset et al. 2010). A combina- 
tion of radiative cooling and three mechanical coolers produces 

^ http : //hea . iki . rssi . ru/rttl58/en/ 
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Table 1. Original sample of 26 clusters. References for cluster information are: (1) Crawford et al. (1995); (2) Appenzeller et al. 
(1998); (3) Ebeling et al. (1998); (4) Bohringer et al. (2000); (5) Ebeling et al. (2002); (6) Ebeling et al. (2007); (7) Kocevski et al. 
(2007); (8) Planck Collaboration VIII (2011); and (9) Planck RTT follow-up programme; see Section 2. The meaning of the notes 
in the right-most column are: (i) cluster has SZ detection smaller than 3cr before source subtraction and smaller than 5cr after source 
subtraction; (ii-a) cluster has a total integrated source flux density greater than 5 mJy within a radius of 3'; (ii-b) cluster has a total 
integrated source flux density greater than 10 mJy within a radius of 10'. "in sample" indicates that the cluster is included in the 
1 1 -cluster sample analysed in this paper 



ESZ cluster name 


Alternative cluster name 


Right Ascension 


Declination 


RA and Dec 


Redshift 


Redshift 


Rejection 






(J2000) 


(J2000) 


reference 




reference 


reason 


PLCKESZG139. 59+24.18 




06:21:58.00 


+74:42:15.7 


8 


0.27 


9 


in sample 


PLCKESZ G167. 65+17.64 


ZwCl 0634.1+4750 


06:38:04.00 


+47:47:53.0 


1,7 


0.17 


1,7 




PLCKESZ G056.80+36.30 


MACS J0717.5+3745 


07:17:30.93 


+37:45:29.7 


6 


0.55 


6 


ii-b 


PLCKESZ G157.43+30.33 


RXJ0748+5941 


07:48:45.60 


+59:41:41.0 


2 


0.55 


2 


in sample 


PLCKESZ G149.70+34.70 


A665 


08:30:58.50 


+65:51:03.7 


3 


0.18 


3 


ii-b 


PLCKESZ G186.40+37.30 


A697 


08:42:57.56 


+36:21:59.3 


3 


0.28 


3 


ii-b 


PLCKESZ G166.13+43. 39 


A773 


09:17:52.97 


+51:43:55.5 


3 


0.22 


3 


in sample 


PLCKESZ G 1 95 . 60+44. 1 


A781 


09:20:26.08 


+30:30:54.0 


3 


0.30 


3 


ii-a 


PLCKESZ G163.70+53.50 


A980 


10:22:28.10 


+50:07:15.6 


3 


0.16 


3 


ii-b 


PLCKESZ G165.08+54. 11 


A990 


10:23:41.83 


+49:08:38.2 


3 


0.14 


3 


in sample 


PLCKESZ G228. 15+75. 19 


MACS J 1149+2223 


11:49:34.30 


+22:23:42.5 


6 


0.55 


6 


in sample 


PLCKESZ G226.24+76.76 


A1413 


11:55:18.24 


+23:24:28.6 


3 


0.14 


3 


in sample 


PLCKESZ G180.60+76.70 


A 1423 


11:57:22.10 


+33:37:55.2 


3 


0.21 


3 




PLCKESZ G229. 60+78.00 


A 1443 


12:01:17.00 


+23:06:18.0 


4 


0.27 


4 


ii-a,b 


PLCKESZ G125.70+53.90 


A1576 


12:36:58.96 


+63:11:26.5 


3 


0.30 


3 


ii-a 


PLCKESZ G121. 11+57.01 




12:59:23.80 


+60:05:24.8 


8 


0.34 


9 


in sample 


PLCKESZ Gl 18.40+39.30 


RXCJ1354.6+7715 


13:54:37.80 


+77:15:34.6 


4 


0.40 


4 




PLCKESZ G067. 23+67.46 


A1914 


14:26:02.15 


+37:50:05.8 


3 


0.17 


3 


in sample 


PLCKESZ G053. 52+59.54 


A2034 


15:10:10.80 


+33:30:21.6 


3 


0.11 


3 


in sample 


PLCKESZ G044.20+48.70 


A2142 


15:58:22.10 


+27:13:58.8 


3 


0.09 


3 


ii-b 


PLCKESZ G097.73+38. 1 1 


A2218 


16:35:52.80 


+66:12:50.0 


3 


0.17 


3 


in sample 


PLCKESZ G086.50+ 15.30 


A2219 


16:40:22.56 


+46:42:32.4 


3 


0.23 


3 


ii-a,b 


PLCKESZ G056.80+36.30 


A2244 


17:02:42.87 


+34:03:42.8 


3 


0.10 


3 


ii-a 


PLCKESZ G055. 60+3 1.90 


A2261 


17:22:27.09 


+32:08:01.7 


3 


0.22 


3 


ii-a,b 


PLCKESZ G086.50+ 15.30 


CIZA J1938.3+5409 


19:38:18.60 


+54:09:33.0 


5 


0.26 


5 




PLCKESZ G077.90- 26.64 


A2409 


22:00:53.03 


+20:57:38.3 


3 


0.15 


3 


in sample 



the temperatures needed for the detectors and optics (Planck 
Collaboration II 2011). Two Data Processing Centers (DPCs) 
check and calibrate the data and make maps of the sky (Planck 
HFI Core Team 2011b; Zacchei et al. 2011). Planck's sensitiv- 
ity, angular resolution, and frequency coverage make it a pow- 
erful instrument for Galactic and extragalactic astrophysics as 
well as cosmology. Early astrophysics results are given in Planck 
Collaboration VIII-XXVI 2011, based on data taken between 
13 August 2009 and 7 June 2010. Intermediate astrophysics re- 
sults are now being presented in a series of papers based on data 
taken between 13 August 2009 and 27 November 2010. 

We note that the Planck maps used for the analysis in this 
paper are not the same as those used in the Early Planck results 
papers. In particular, we stress that both the data and the analysis 
techniques employed for this study are not the same as was used 
to construct the Planck ESZ catalogue (Planck Collaboration 
VIII 2011). A later version of the data has been used for the 
analysis here and we use different analysis techniques for rea- 
sons which will be explained in Section 5.1. However, as part of 
our suite of internal tests, we have repeated our analysis on the 
older version of the Planck data which was used to derive the 
ESZ catalogue and we find excellent agreement. 

In Fig. 1 we present maps of the dimensionless Compton-y 
parameter for each of the clusters as estimated from the Planck 



data. The y-parameter is related to the observed brightness as a 
function of frequency (v) by 

^r^(v)=yg(v), (1) 

^ CMB 

where AT^sz is the brightness fluctuation due to the thermal 
SZ (tSZ) effect and Tcmb is the temperature of the CMB, 
which we take to be 2.7255 + 0.0006 (Fixsen 2009). The 
function, g{v) is the frequency dependence of the SZ effect 
(Sunyaev & Zeldovich 1972). For the HFI channel frequencies 
(100,143,217,353,545 and 857 GHz), g{v) takes the values 
(-4.03,-2.78,0.19,6.19, 14.47 and 26.36 KcuB/y)- The maps 
of Fig. 1 were estimated from HFI channel data taken be- 
tween 13 August 2009 and 27 November 2010, corresponding 
to slightly more than 2.5 full-sky scans. The measured noise lev- 
els on the Planck frequency channel maps are listed in Table 2 
(Planck HFI Core Team, in prep.). The tSZ signal reconstruc- 
tions were performed using the MILCA method (Hurler et al. 
2010 and references therein) on the six Planck all-sky maps 
from 100 GHz to 857 GHz. MILCA (Modified Internal Linear 
Combination Algorithm) is a component separation approach 
aimed at extracting a chosen component (here the tSZ signal) 
from a multi-channel set of input maps. It is based on the well 
known ILC approach (e.g., Eriksen et al. 2004), which searches 
for the linear combination of the input maps that minimises the 
variance of the final reconstructed map, while imposing spec- 
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Table 2. Noise levels per A^side = 2048 pixel for each Planck 
frequency band. 



Freq. (GHz) 


100 


143 


217 


353 


545 


857 


Noise (//K) 


80 


34 


56 


167 


1414 


170200 



tral constraints. For our cluster SZ reconstructions, we applied 
MILCA using two constraints, the first one to preserve the tSZ 
signal and the second one to remove CMB contamination in the 
final tSZ y-map. In addition, to compute the weights of the lin- 
ear combination, we have used the extra degrees of freedom to 
minimise residuals from other components (two degrees) and 
from the noise (two degrees). The noise covariance matrix was 
estimated from jack-knife maps. The final y-maps have an ef- 
fective resolution of 10 arcmin. Note that, in general, the prop- 
erties of the foreground emission depend on both the position 
on the sky and on the frequency of observation. We have there- 
fore allowed the weights to vary as a function of both position 
and frequency. We have confirmed using simulations that such 
an approach maximises the signal-to-noise and minimises the 
bias in the extraction of the tSZ signal. We emphasise that the 
MILCA SZ reconstructions presented in Fig. 1 are intended for 
visual examination and qualitative assessment of the cluster sig- 
nals only. Our quantitative analysis of the Planck data, which 
we use to compare with the AMI results, is based mainly on 
the PowellSnakes (Carvalho et al. 201 1) SZ extraction algorithm 
(see Section 5). 

4. Description of AlVII data 

AMI comprises two arrays, the Small Array (SA) and the 
Large Array (LA), located at the Mullard Radio Astronomy 
Observatory near Cambridge. The SA consists of ten 3.7-m di- 
ameter equatorially-mounted antennas, with a baseline range 
of ^ 5-20 m and synthesised beam (resolution) of around 3'. 
The LA consists of eight 13-m diameter antennas with a base- 
line range of ^ 20-100 m and synthesised beam of around 30". 
Both arrays observe Stokes I + Q '\n the band 13-18 GHz, each 
with system temperatures of about 25 K. Note that AMI defines 
Stokes Q and U with respect to celestial north. The backends are 
analogue Fourier transform spectrometers, from which the com- 
plex signals in each of eight channels of 750-MHz bandwidth 
are synthesised, and the signals in the adjacent channels are cor- 
related at the ^10 % level. Further details of the instrument are 
given in AMI Consortium: Zwart et al. (2008). 

SA pointed observations of our cluster sample were taken 
during 2007-2011. The observation lengths per cluster before 
any flagging of the data are presented in Table 4; the noise levels 
on the SA maps reflect the actual observation time used. The SA 
observations were made with single pointings interspersed with 
a phase calibration source, while the LA observations were made 
in a 6 1 + 1 9-point raster mode configuration with 4' spacing . This 
consisted of 61 pointings arranged in a hexagonal grid, with grid 
points separated by 4' with further observations of the central 19 
pointings designed to increase the sensitivity at the centre of the 
field. In this mode the integration time on the area < 7.5' from 
the cluster centre is twice as long as the area > 7.5' away, so as 
to provide a better match to the primary beam sensitivity of the 
AMI SA. Phase calibrators were chosen from the Jodrell Bank 
VLA Astrometric Survey (JVAS, Patnaik et al. 1992) on the ba- 
sis of proximity (< 2° for the AMI LA, < 8° for the AMI SA) 



Table 3. Assumed I + Q flux densities of 3C286 and 3C48 over 
the commonly-used AMI band, and the full width at half maxi- 
mum of the LA primary beam (approximate field of view, ©la) 
for each channel. 



Channel 


y/GHz 


53C286/jy 


53«8/Jy 


©la/ arcmin 


3 


13.9 


3.74 


1.89 


6.08 


4 


14.6 


3.60 


1.78 


5.89 


5 


15.3 


3.47 


1.68 


5.70 


6 


16.1 


3.35 


1.60 


5.53 


7 


16.9 


3.24 


1.52 


5.39 


8 


17.6 


3.14 


1.45 


5.25 



and 15 GHz flux density (> 0.2 Jy for the AMI LA, > 0.7 Jy for 
the AMI SA). The JVAS is based on observations made with the 
VLA in "A" configuration (Condon & Broderick 1985, 1986; 
Condon et al. 1989; White & Becker 1992). 

The reduction of the AMI data was performed using a 
dedicated software tool reduce. This is used to apply path- 
compensator and path-delay corrections, to flag interference, 
shadowing and hardware errors, to apply phase and amplitude 
calibrations and to Fourier transform the correlator data readout 
to synthesise the frequency channels, before outputting to disk. 

Flux calibration was performed using short observations of 
3C48 and 3C286 near the beginning and end of each run. The 
assumed I+Q flux densities for these sources in the AMI chan- 
nels are listed in Table 3 and are consistent with Baars et al. 
(1977). As Baars et al. measure / and AMI measures I+Q, 
these flux densities include corrections for the polarisation of 
the sources. An amplitude correction is also made for the inter- 
vening air mass during the observation. Flux calibration is ex- 
pected to be accurate to ^3% for the AMI SA and ^5% for 
the AMI LA. After phase calibration, the phase of both arrays 
over one hour is generally stable to 5° for channels 4-7, and to 
10° for channels 3 and 8. (Channels 1 and 2 are generally not 
used for science analysis as they tend to suffer from interference 
problems.) 

Maps were made using the Astronomical Image Processing 
System (aips, Greisen 2003) from each channel of the AMI SA 
and LA; however here we present only the combined-channel 
maps of the SA and LA observations. The aips task imean was 
used on the LA individual maps to attach the map noise to the 
map header imean fits a Gaussian to the histogram of the map 
pixels (ignoring extreme pixels that might be due to sources) and 
uses the standard deviation of the fitted Gaussian as a measure 
of the random noise in the data. The aips task flatn was then 
used to form a mosaiced image from the multiple pointings. Data 
from the pointings were primary beam corrected using parame- 
ters listed in Table 3 and weighted accordingly when combined. 
The AMI SA combined-channel map noise and the LA map 
noise are given in Table 5. The raw uv data for all good obser- 
vations were concatenated together to make a visibility data file 
for each channel. All maps were made using natural uv weight- 
ing and all images were cleaned to three times the thermal noise 
with a single clean box encompassing the entire map. The data 
were also binned into bins of width 40/1. This reduced the size 
of the data to a manageable level without adversely affecting the 
subsequent inference of cluster properties. Fig. 2 shows a typ- 
ical example of the SA synthesised beam, in this case for the 
observations of A2218. 
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Fig. 1. Maps of the dimensionless Compton-y parameter (equation 1) as estimated from the Planck observations using the MILCA 
algorithm. The maps have an effective resolution of 10 arcmin. The clusters are ordered in terms of increasing redshift, from top left 
to bottom right. Each panel shows a 100' x 100' region. 
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Fig. 2. SA synthesised beam for the A2218 observations. 
Contours start at 6 % and increase linearly by 3 % per contour 
Contours drawn as red dashed lines are negative. The synthe- 
sised beams for the other cluster observations are qualitatively 
similar. 



As contamination from radio sources at 16 GHz tends to be 
significant, removing or modelling this emission accurately can 
often be essential to recover SZ decrements from AMI maps. To 
address this issue, we use the AMI-developed source extraction 
software, sourcefind (Waldram et al. 2003; AMI Consortium: 
Franzen et al. 2011) to determine the position, flux density and 
spectral index of the radio sources with flux density > 3.5crLA 
on the CLEANED LA continuum maps, where ctla is the LA ther- 
mal noise. Spectral indices were fit with a Markov Chain Monte 
Carlo (MCMC) method using LA maps for all six channels and 
assuming that source flux densities follow a power-law relation 
of 5 oc v^" for the AMI frequencies. These source parameter es- 
timates are subsequently used as priors in our Bayesian analysis 
of the cluster SZ signals (Section 5.2). Tables 4 and 5 summarise 
the observational details of our cluster sample and Figs. 3 and 4 
present the maps of the AMI observations of these clusters be- 
fore and after source subtraction, respectively. Once again, as 
with the reconstructed Planck maps presented in the previous 
section, the AMI maps presented in Figs. 3 and 4 are intended 
for visual examination and qualitative assessment of the cluster 
signals only. Our quantitative analysis of the AMI data is de- 
scribed later in Section 5.2. 



5. Analysing the SZ signal 

The SZ surface brightness (the Compton-y parameter, cf., equa- 
tion 1) is proportional to the line of sight integral of the electron 
pressure. 



y = 



f 

•J —C 



P.ir)dl, 



(2) 



where Pe{r) is the electron pressure at radius r, crj is the 
Thomson scattering cross-section, nif. is the electron mass, c 
is the speed of light and dl is the line element along the line 
of sight. In this context, Nagai et al. (2007) analysed the pres- 
sure profiles of a series of simulated clusters (Kravtsov et al. 



2005) as well as a sample of relaxed real clusters presented in 
Vikhlinin et al. (2005, 2006). They found that the pressure pro- 
files of all of these clusters could be described by a generali- 
sation of the Navarro, Frenk, and White (NFW, Navarro et al. 
1997) model used to describe the dark matter halos of simulated 
clusters. Assuming spherical geometry, the GNFW pressure pro- 
file (Nagai et al. 2007) reads 



Pe(r)^Po\- 



(y-l3)la 



(3) 



where Pq is the overall normalisation coefficient of the pres- 
sure profile and is the scale radius. It is common to define 
the latter in terms of rsoo, the radius at which the mean den- 
sity is 500 times the critical density at the cluster redshift, and 
to define the gas concentration parameter, C500 = '"50o/'"s- The 
parameters {a,/3,y) describe the slopes of the pressure profile 
at r ^ Ts, r > r^, and r <s: respectively. In order to retain 
consistency between the Planck and AMI analysis pipelines, we 
follow Arnaud et al. (2010) (see also Planck Collaboration VIII 
2011; Planck Collaboration X 2011) and fix the values of the 
gas concentration parameter and the slopes to be (cjoo, q:,/3, 7) = 
(1.156, 1.0620,5.4807,0.3292). These values describe the "uni- 
versal pressure profile", derived from XMM-Newton observa- 
tions of the REXCESS cluster sample (Bohringer et al. 2007), 
and from three different sets of detailed numerical simulations 
by Borgani et al. (2004), Piff^aretti & Valdarnini (2008), and 
Nagai et al. (2007), which take into account radiative cooling, 
star formation, and energy feedback from supernova explosions. 
In Section 7, we will relax these restrictions for a subset of our 
cluster sample and will include information from X-ray obser- 
vations of individual clusters in our analysis. We note that the 
profile of equation (3) has recently been used to analyse SZ data 
from the South Pole Telescope (SPT, Plagge et al. 2010) in addi- 
tion to the Planck survey data (Planck Collaboration VIII 201 1; 
Planck Coflaboration X 201 1). 

The integral of the y parameter over the solid angle Q sub- 
tended by the cluster is denoted by Y^z, and is proportional to 
the volume integral of the gas pressure. It is thus a good indi- 
cator of the total thermal energy content of the cluster and its 
mass (e.g., Bartlett & Silk 1994). The determination of the nor- 
malisation and the slope of the Fsz - M relation has therefore 
been a major goal of studies of the SZ effect (da Silva et al. 
2004; Nagai 2006; Kravtsov et al. 2006; Plagge et al. 2010; 
Amaud et al. 2010; Andersson et al. 201 1 ; Planck Collaboration 
VIII 201 1; Planck Collaboration IX 201 1; Planck Coflaboration 
X 2011; Planck Coflaboration XI 2011; Planck Coflaboration 
XII 201 1). In particular, Andersson et al. (201 1) investigated the 
YsL - Yx scaling relation within a sample of 15 clusters observed 
by SPT, Chandra and XMM-Newton and found a slope of close 
to unity (0.96 + 0.18). Similar studies were carried out by Planck 
Collaboration XI (2011) using a sample of 62 nearby (z < 0.5) 
clusters observed by both Planck and by XMM-Newton. The re- 
sults are consistent with predictions from X-ray studies (Arnaud 
et al. 2010; Andersson et al. 2011). These studies at low red- 
shifts, where data are available from both X-ray and SZ obser- 
vations of galaxy clusters, are crucial to calibrate the Ysz - M 
relation, as such a relation can then be scaled and used to deter- 
mine masses of SZ selected clusters at high redshifts in order to 
constrain cosmology. 

The integrated y parameter (Fsz) adopting a spherical geom- 
etry Ksph, is given by 



CTj 

nif^c' 



Jo 



P^(,r')Anrdr'. 



(4) 
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Fig. 3. AMI maps before source subtraction. The clusters are ordered as in Fig. 1, in terms of increasing redshift. Black solid lines 
represent positive contours and red dashed lines indicate negative contours. The contours increase linearly from +2crsA to +10crsA 
where ctsa is listed in Table 5 for each cluster Each map covers a region approximately 30' x 30' and the resolution is around 3'. 

7 



Planck and AMI Collaborations: Comparison of Planck and AMI Sunyaev-Zeldovich measurements for 1 1 galaxy clusters 









' '^j 1 u ' ■ 


o 


o 




- % 




^ C'-, - 


Q 








o 




r 




. , ^ ■ 



53,7 63.6 53.5 53,4 
Gal, longitude (degs) 



226,4 226,3 226,2 226.1 226 
Gal. longitude (degs) 



165.2 165.1 165 164. £ 
Gal. longitude (degs) 




78 77.9 77.8 77.7 
Gal. longitude (degs) 





67.4 67.3 67.2 87.1 
Gal. longitude (degs) 



O 



o 



97.9 97.8 97.7 97.8 
Gal. longitude (degs) 



PLCKESZ G139. 59 + 24. IB 



PLCKESZ G121. 11+57.01 











§ 




C 






■ 















166,3 166,2 166,1 166 165,9 
Gal, longitude (degs) 



^ (~l 



'"^ 







139.8 139.7 139.6 139.5 
Gal. longitude (degs) 




121.3 121.2 121.1 121 120.9 
Gal. longitude (degs) 



MACS J1149+2223 



RXJ0748+5941 



1 S 




226.4 228,3 228,2 226.1 228 
Gal, longitude (degs) 




" 



157.6 157.5 157.4 157.3 
Gal. longitude (degs) 



Fig. 4. AMI maps after subtraction of radio point sources. The solid black lines represent positive contours and the dashed red lines 
indicate negative contours. The cluster ordering, contour levels and resolution are the same as in Fig, 3, 
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Table 4. Details of the AMI observations of our cluster sample. 



Cluster 


Observing dates 


AMI SA observing 


LA observing 


AMI SA phase 


LA phase 






time (hours) 


time (hours) 


calibrator 


calibrator 


A2034 


9010 Tan Mar Dpr 
aaj L\j jaii, iviai, Ljc^ 


38 

JO 


93 


J i JUUT J 1 jyj 


Tl ^044-3940 




9007 Mar NTnir 9010 Mar 
I iVldl, i>UV, Z.UiU iVldl 


40 


10 




Tl 1 S04-941 7 




9007 Mr»\7- 900Q Ar»r 


J 1 


17 
1 / 




Tl 01 S4-ZLQ9/^ 


A2409 


2007 Mar, Apr, May; 2009 Dec 


38 


21 


J2225+2I18 


J2200+2137 


A1914 


2008 May; 2009 Jan, Jun 


32 


35 


J1419+382I 


JI419+3821 


A2218 


2008 Jan; 2009 Feb, Jun 


30 


12 


J 1642+6856 


JI623+6624 


A773 


2007 Sep, Oct; 2009 May 


40 


19 


J0903+465I 


J0929+5013 


MACS Jl 149+2223 


2010 Apr, May, Nov, Dec 


38 


11 


Jl 150+2417 


Jl 150+2417 


RXJ0748+5941 


2011 Feb 


45 


35 


J0753+5352 


J0737+5941 


PLCKESZG139.59+24.18 


2010 Nov 


52 


37 


J0639+7324 


J0639+7324 


PLCKESZGI21. 11+57.01 


2011 Jan 


40 


67 


J 1302+5748 


J1302+5748 



Table 5. Details of radio point sources and the thermal noise levels for the AMI observations of the cluster sample. Here ctsa and 
(Tla refer to the thermal noise levels reached in the LA and SA maps, respectively. The S/N values are calculated by dividing the 
peak flux values by cr^p^. 



Cluster 


Number of LA 
3.5 (Tla sources 


total source flux 
on the map (mJy beam"') 


CTsA 

(mJy beam"') 


ctla 
(mJy beam"') 


S/N 


A2034 


13 


32.5 


0.08 


0.08 


10 


A1413 


20 


47.6 


0.09 


0.12 


12 


A990 


20 


56.3 


0.11 


0.10 


8 


A2409 


14 


49.2 


0.11 


0.10 


11 


A1914 


14 


31.5 


0.12 


0.10 


11 


A2218 


9 


31.5 


0.11 


0.15 


12 


A773 


6 


12.8 


0.12 


0.12 


6 


MACS J 1149+2223 


13 


44.6 


0.09 


0.14 


17 


RXJ0748+5941 


15 


17.1 


0.07 


0.06 


12 


PLCKESZ G 1 39.59+24. 1 8 


13 


13.4 


0.09 


0.06 


11 


PLCKESZG121. 11+57.01 


19 


25.4 


0.07 


0.06 


10 



Following Arnaud et al. (2010), we consider the radius of Sr^oo 
as the boundary of the cluster where the pressure profile flattens, 
and we use this boundary to define the total volume integrated 
SZ signal, Ftot- 

In the simplest case, where a, /3, y, and C500 in equation (3) 
have fixed values, our cluster model depends only on four pa- 
rameters: Xc and which define the projected cluster position 
on the sky and Pq and in the pressure profile (equation 3). 
In this paper, we define clusters in terms of the parameter set 
©c = {xc,yc, 0s = rJDA, Ye = Ytot/Dj), where Da is the 
angular-diameter distance to the cluster (Planck Collaboration 
VIII 201 1), and we determine the model parameter Pq by evalu- 
ating equation (4) at r = Srsoo. To calculate Da we assume a flat 
Universe with matter density f2„, - 0.27 and Hubble constant 
Ho = 70.4kms"iMpc"' (Komatsu et al. 2011). 

We adopt an exponential prior for 6*, and a power-law prior 
for Yg to analyse both the Planck and AMI data (Carvalho et al. 
2011). The prior on 0^ is ^e"'"'^ for 1.3' < 0^ < 45' and zero 
outside this range, with A = 0.2. The prior on Yg is Yg" for 

5.0 X lO-^arcmin^ < Yg < 0.2 arcmin" and zero outside this 
range, with a - 1.6. These priors have already been used in 
Planck detection and extraction algorithms to identify and char- 
acterise compact objects buried in a diffuse background (Planck 
Collaboration VIII 2011). For the cluster position, however, in 
order to ensure that we are comparing integrated SZ fluxes cen- 
tred on identical positions on the sky, we performed an initial 



analysis of the AMI data using a Gaussian prior centred on the 
cluster phase centre and with a standard deviation of 1 ' in order 
to find the best-fitting cluster coordinates. We then fixed the clus- 
ter position to these best-fitting coordinates in the subsequent 
analysis of the Planck data and also in a subsequent re-analysis 
of the AMI data. 

5.1. Analysis of Planck data 

The analysis of the Planck data was performed using 
PowellSnakes (PwS), which is a Bayesian package for discrete 
object detection, photometry and astrometry, as described in 
Carvalho et al. (2009, 2011). PwS is part of the Planck HFI 
pipeline and is regularly used to produce catalogues of objects 
(Planck Collaboration VII 201 1) and to measure and characterise 
the SZ signal (Planck Coflaboration VIII 2011). Note that we 
have chosen to use PwS as our primary SZ extraction algorithm 
for the Planck analysis in this study as PwS naturally returns the 
full posterior distribution in the Y^oq - 6500 2D parameter space. 
It is thus naturally suited for combining with the AMI results 
to produce joint constraints. We will also present a comparison 
with results obtained using the Matched Multi-Filter algorithm 
(MMF3; Melin et al. 2006) which was the reference algorithm 
adopted for the production of the Planck ESZ catalogue (Planck 
Collaboration VIII 2011). However, in its current implementa- 
tion, the MMF3 algorithm does not produce any information on 
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the correlation between the two cluster parameters of interest 
(Fsoo and 6500) and so producing joint constraints as obtained 
from AMI and Planck via MMF3 is currently not possible. 

The analysis of the Planck data using the PwS algorithm pro- 
ceeded as follows. For each cluster, flat patches (14.7° x 14.7°; 
512x512 pixels) were created using a gnomonic projection, cen- 
tred on the targeted cluster, for each of the Planck HFI chan- 
nels. By operating on such a large patch of sky enough statis- 
tics are collected in order to produce a smooth cross-channel co- 
variance matrix. The position of each cluster was assumed to be 
known precisely (we adopted delta-function priors at the AMl- 
determined position) as described in Section 5. 

The data model for a single isolated cluster located in the 
centre of the patch is then described by 



d(x) = }'/r(0s,x) + n(x). 



where x is sky position, d(x) is a vector containing the data, n(x) 
is the background composed of instrumental noise plus all other 
astronomical components except the SZ signal, / is a vector con- 
taining the SZ surface brightness at each frequency, Y is the total 
integrated Comptonisation parameter, 0s is a parameter control- 
ling the cluster radial scale and r(0s, x) is the convolution of the 
canonical GNFW model integrated along the line of sight with 
the Planck beam at that channel. It is assumed that the back- 
ground is a realisation of a stationary Gaussian random field. 

A direct computation of the likelihood is very expensive. 
Therefore, PwS instead computes the likelihood ratio of two 
competing models describing the data: a cluster is present {H\); 
or no cluster is present {Hq). Note that the latter hypothesis does 
not contain any parameters and therefore only multiplies the tar- 
get likelihood in Hi by a constant. The representation of the like- 
lihood ratio in real space reads 



In 



-Cum 



i'r-'[;P/i/)7(-J/;0s)] 
iF2 2Q,7(i7)PF(i7;0,)| 



x=0 

2 



(6) 



where rf is the spatial frequency (the conjugate variable to x) 
and ;F ' [. . .]x denotes the inverse Fourier transform of the quan- 
tity in brackets, evaluated at the point x. We have also de- 
fined the quantities 'Pjirf) = d'{i])N \ri)ipii]) and Qijirj) s 
\lf'-(ri)N^^{ri)ip jirf), in which the vector i^,(j/) has the compo- 
nents (i/f,)v = Bv(i/)(f,)v, with V labeling frequency channels and 
Byirf) is the beam transfer function. The quantity t(-j/; Q^) is the 
Fourier transform of V{Q^,x) and the matrix N{Tf) contains the 
generalised noise cross-power spectra. We refer the interested 
reader to Carvalho et al. (2009, 2011) for further technical de- 
tails on the PwS algorithm. 

The cross-channel covariance matrix is computed directly 
from the pixel data, by averaging the Fourier modes in concen- 
tric annuli. This operation is only possible because of the as- 
sumed isotropy of the background. To reduce the contamination 
of the background by the signal itself, the estimation of the co- 
variance matrix is performed iteratively. After an initial estimate, 
all detected clusters in the patch are subtracted from the data us- 
ing their best fit values and the covariance matrix is re-estimated. 
To enforce our assumption of a single source in the centre of the 
patch, PwS removes from the data all other detections with SNRs 
higher than our target cluster to reduce possible contamination 
of the signal by power leakage from nearby objects. Bright point 
sources are masked or subtracted from the maps as part of a pre- 
processing routine run prior to the production of the flat patches. 



To construct the joint posterior distributions of {Y,0s), we 
have used the set of priors as described in Section 5. To draw the 
posterior distribution manifold, PwS grids the parameter space 
using a uniformly spaced lattice of (256 x 256) cells, appropri- 
ately chosen to enclose all posterior regions significantly differ- 
ent from zero. 

Since the LFI channels of Planck have relatively coarse res- 
olution, the use of LFI bands in current implementations of the 
extraction algorithms results in beam dilution of the SZ signal 
and thus decreases the S/N for the detected clusters (Planck 
Collaboration VIII 2011). This can potentially be improved in 
the future with modifications to the algorithms but for the pur- 
poses of the present study, we use only the HFI data. 



(5) 5.2. Analysis of AMI data 



An interferometer like AMI operating at a frequency v measures 
samples from the complex visibility plane 7v(u). These are given 
by a weighted Fourier transform of the surface brightness /v(x), 
namely 



/v(u) 



Ay(x)Iy(x) exp(27r/u ■ x)dx. 



(7) 



where x is the position relative to the phase centre, Ay(x) is the 
(power) primary beam of the antennas at observing frequency v 
(normalised to unity at its peak) and u is the baseline vector in 
units of wavelength. In our model, the measured visibilities are 
defined as 



Vy(u)^Sy{u)+Ny{u), 



(8) 



where the signal component, 5'v(u), contains the contributions 
from the SZ cluster and identified radio point sources, whereas 
the generalised noise part, Ny{u), contains contributions from a 
background of unsubtracted radio point sources, primary CMB 
anisotropics and instrumental noise. 

We assume a Gaussian distribution for the generalised noise. 
This then defines the likelihood function for the data 



£(©) = ^exp(-i;^2 
Zn \ / 



(9) 



where is the standard statistic quantifying the misfit between 
the observed data D and the predicted data £)''(©), 



= Y^iDy - DPy/(Cyyr\Dy, - D^,), 



(10) 



where y and V are channel frequencies. Here C is the generalised 
noise covariance matrix 



C_ ^rec , ^CMB , ^cc 
— *-'r,v' v.v' V. 



(11) 



The first term on the right hand side of equation (11) is a di- 
agonal matrix with elements cr^,jSijSyy', where cry, is the rms 
Johnson (receiver) noise on the ;th element of the data vector Dy 
at frequency v. The second term denotes the noise due to pri- 
mordial CMB anisotropics and contains significant off-diagonal 
elements, both between visibility positions and between fre- 
quencies. This term can be calculated from a given primary 
CMB power spectrum C™®(v) following Hobson & Maisinger 
(2002); note that in intensity units the CMB power spectrum is a 
function of frequency. To calculate this term, we adopt the best- 
fitting CMB power spectrum to the WMAP 7-year data (Komatsu 
et al. 2011). The third term on the right hand side of equa- 
tion (11) is the source confusion noise, which accounts for re- 
maining unresolved radio sources with flux densities less than 
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the flux limit (S'lim) of the AMI observations and which remain 
after high resolution observation and subtraction. We estimate 
this term assuming that sources are randomly distributed on the 
sky, in which case we can describe the source confusion noise 
with a power spectrum calculated as 

C™"f(v) = f S\y{S)AS, (12) 
Jo 

where ny{S) - dNy{> S)/dS is the differential source count at 
frequency v as a function of flux density S . We use the source 
counts as measured by the IOC survey (AMI Consortium: Davies 
et al. 201 1) for our calculation. The limiting flux density for the 
integration (5 Mm) is determined from the noise in the LA maps 
and is different for each cluster, but is typically in the range 0.2- 
0.5 mJy. 

The normalisation factor Zn in equation (9) is given by 
Zn = {27:f''-''>'^\C\"\ (13) 

where Nyis is the total number of visibilities. Further details on 
our Bayesian methodology, generalised noise model, likelihood 
function and resolved radio point-source models are given in 
Feroz & Hobson (2008) and Feroz et al. (2009a,b). 

Radio sources detected in the LA maps were modeled by 
four source parameters, 0s = (xs , , S o , a), where jCj and 
refer to the right ascension and declination of radio sources, re- 
spectively, while S and a are the flux density and spectral index 
of the radio source at the central frequency, vq. As mentioned in 
Section 4, this modelling is necessary because of source variabil- 
ity and some difficulty with inter-an^ay calibration. Therefore, 
the properties of point sources detected at > S.Sctla by the 
LA were used as priors when modelling the SA data. We used 
a delta-function prior on the position of the source since the 
resolution of the LA is around three times that of the SA. We 
used Gaussian priors on the source flux densities, with LA (in- 
tegrated, where applicable) flux densities generating the peak of 
the prior, and the Gaussian crs were set to 40 % of the source flux 
densities. Spectral index (a) priors were also set as Gaussians, 
with cr equal to the error on the spectral index fit. This is be- 
cause for sources with high signal-to-noise ratio, the determina- 
tion of the spectral index is dominated by the AMI frequency 
channel mean and the error on a is Gaussian. For sources with 
low signal-to-noise ratio which just meet our continuum detec- 
tion threshold, the spectral index probability distribution is dom- 
inated by the prior, which is determined from the IOC survey 
(AMI Consortium: Davies et al. 201 1). 

5.3. Results 

Fig. 5 presents the 2D marginalised posterior distributions in the 
^500 - ^500 plane and Table 6 summarises the mean and the dis- 
persion of these two parameters for each cluster, as estimated 
from the Planck and AMI data, respectively. Note that in Fig. 5 
the inner and outer contours show the areas enclosing 68 % and 
95 % of the probability distributions. Estimates of 6500 as derived 
from X-ray observations, and which were included in the Planck 
ESZ catalogue (Planck Collaboration VIII 2011), are also indi- 
cated in the figures for comparison. 

Recall that in this figure the clusters are ordered in terms of 
increasing redshift. The constraints from both Planck and AMI 
demonstrate a strong cluster size-integrated Compton parameter 
(0500 - i'soo) degeneracy/correlation in all cases. It is important to 
account for such effects when attempting to use the SZ signal to 
estimate cluster masses (da Silva et al. 2004; Arnaud et al. 2007). 



We also note that the Planck constraints appear to be weaker for 
high redshift clusters, which can generally be understood as a 
resolution effect - Planck's relatively poor resolution (e.g., as 
compared with AMI) means it has difficulty resolving and thus 
estimating the parameters of clusters with small angular extent - 
and high-redshift clusters are likely to be smaller in angular size. 
ami's increased resolution, on the other hand, means that it can 
still constrain the sizes of these high-redshift, small-angular size 
clusters. 

For three clusters (A1413, A1914, and PLCKESZ 
G139.59h-24.18), the AMI and Planck constraints are clearly 
discrepant. On the other hand there is significant overlap in 
the posterior distributions for the remaining eight clusters. 
However, taking our cluster sample as an ensemble, there is 
some evidence that the cluster parameter estimates derived from 
the AMI data are systematically lower than those derived from 
the Planck data (i.e., AMI is finding the clusters to be fainter 
and smaller in angular extent compared to what the Planck data 
indicate). 

In addition, the F500 - ^^500 degeneracies are significantly dif- 
ferent for the Planck and AMI constraints, with the degeneracies 
of the AMI constraints being generally steeper than the Planck 
ones. This arises because of the interplay between the angular 
size-redshift relation and the differing angular scales that AMI 
and Planck are sensitive to, as well as the very different observa- 
tional techniques and frequencies used by the two instruments. 

In the cases where the Planck and AMI-derived constraints 
are compatible with one another, we also overplot the joint con- 
straints obtained from multiplying the Planck and AMI posteri- 
ors. In many cases, the resulting joint constraints are far tighter 
than either analysis alone which is a direct result of the differ- 
ing parameter degeneracies for the two instruments, as described 
above. The marginalised constraints from this combined analysis 
are also presented in Table 6. 

In Fig. 5, we have also over-plotted the constraints as 
obtained from the Planck data using the MMF3 algorithm. 
Comparing these results with the PwS Planck results, we see 
good agreement in most cases, although there may be a tendency 
for the MMF3 estimates to be systematically brighter and larger 
than the PwS results. However, it is clear that our broad con- 
clusions regarding the general levels of agreement between the 
Planck and AMI results remain unchanged if we consider the 
MMF3 results in place of the PwS constraints. 

In Fig. 6, we plot the Planck determined integrated 
Compton- y parameter versus the Compton- 7 parameter as de- 
rived from the AMI data. Note that, for this correlation plot, 
we have fixed the cluster scale size to be that determined from 
X-ray observations (as indicated by the red arrows in Fig. 5). 
(Three of the clusters have no reported X-ray size so only eight 
of the 11 clusters contribute to this correlation analysis.) The 
measured correlation coefficient is 0.79 and the best-fitting lin- 
ear relationship has a slope of 1.18 + 0.07, again indicating that 
the Planck SZ fluxes appear to be systematically larger that the 
AMI derived fluxes. We have also repeated this analysis fixing 
the cluster size to both the f /a«cA:-determined size and the AMI- 
determined size. In both cases we see the same general trend, 
with the Planck Compton- 7 parameter being consistently larger 
than the AMI-derived value. 

In summary, our results suggest a systematic difference be- 
tween the Planck and AMI measurements of the SZ signal com- 
ing from our cluster sample. Such a systematic difference could 
be an indicator of a shortcoming in some part of our analysis 
and could have important implications for performing cosmo- 
logical studies with larger samples of SZ clusters. For example. 
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Fig. 5. Recovered Planck and AMI 2D posterior distributions in the Fgoo - ^500 plane. Blue contour plots are the results from the 
AMI analysis and purple contours show the Planck results (specifically using the PwS method). Red arrows show the values of 
^500 as determined from X-ray measurements of these clusters where available. The grey points with error bars show the MMF3 
Planck results. The inner and outer contours in each set indicate the areas enclosing 68 % and 95 % of the probability distribution, 
while the MMF3 error-bars indicate the 1 cr uncertainties. Where the recovered AMI and PwS Planck constraints are consistent, the 
joint constraints are also indicated by the heavy black contours. In cases where the contours do not close at the lower ends of the 
parameter ranges, the corresponding constraints represent upper limits only. 

12 



Planck and AMI Collaborations: Comparison of Planck and AMI Sunyaev-Zeldovich measurements for 1 1 galaxy clusters 



Table 6. Recovered mean and dispersion values for 6500 and F500 for the 1 1 clusters. Where consistency is found between the Planck 
and AMI measurements, the joint constraints are also given. The cluster redshift and signal-to-noise of the PwS detections are also 
listed. 
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Fig. 6. Comparison of the integrated Compton-F parameters ob- 
tained from the Planck and AMI fits, when the cluster size is 
fixed to that determined from X-ray observations. The one-to- 
one relation is denoted with the dashed line. The best-fitting lin- 
ear relation is plotted as the unbroken line. The slope of this 
latter relation is 1.18 + 0.07 and the correlation coefficient is 
0.79. Note that the same general behaviour (slope > 1) is also 
observed when we fix the cluster size to be that determined from 
either the Planck or the AMI SZ observations. 



the observed systematic could indicate that the way the clusters 
are being modeled in the analysis (e.g., the fixed GNFW pro- 
file adopted and/or the assumption of spherical symmetry) is not 
flexible enough to describe both the Planck and AMI results si- 
multaneously. If this were to be the source of the discrepancy 
then such effects would need to be accounted for in future cos- 
mological studies. However, before considering such an expla- 
nation, it is important to first consider if possible instrumental 
and/or astrophysical systematic effects could be responsible for 
the results we have found. We now turn to simulations to inves- 
tigate the potential impact of such effects. 



6. Simulations 

In order to test the SZ signal extraction techniques used and to 
investigate whether the systematic discrepancy observed in the 
real data is due to unaccounted-for astrophysical foregrounds, 
instrumental systematics or data-analysis induced biases, we 
have conducted detailed simulations of both the Planck and AMI 
experimental setups. 

For each of the 1 1 clusters in our sample, to create an input 
SZ signal for the simulations, we simulated a cluster SZ signal 
using the GNFW pressure profile (equation 3), with input param- 
eters based on the best-fitting Y500 and 6500 values from an anal- 
ysis based on intermediate Planck maps, which are, in practice, 
close to the best- fitting Planck parameters quoted in Table 6. 

6.1. Planck simulations 

For Planck, the main worry in terms of astrophysical systemat- 
ics is probably thermal emission from dust in the Galaxy. The 
Planck simulation ensemble comprised of CMB and noise re- 
alisations and a fixed foreground dust template, produced by 
re-scaling the Planck 857 GHz channel map to the other HFI 
frequencies and reconvolving so as to apply the appropriate 
beam for each channel. The dust template assumed a modified 
blackbody spectrum with emissivity y6 = 1.8 and temperature 
T - 18 K. The beams were assumed to be Gaussian with the 
appropriate mean FWHM for each channel as calculated by the 
FEBeCoP algorithm (Planck HFI Core Team 201 lb). 

The noise component of the simulations was generated us- 
ing the Springtide destriping pipeline (Ashdown et al. 2007). 
This pipeline creates realisations of the noise in the nominal mis- 
sion time-ordered data streams, compresses them to rings and 
destripes the rings to produce noise maps. It is assumed that the 
noise is uncorrelated between rings and that in each ring it is 
drawn from a power spectrum. For these simulations, the noise 
power spectrum used was the mean of the ring -by-ring spectra. 
In turn, these were determined by applying the noise estimation 
pipeline (Planck HFI Core Team 201 lb) to the exact same ver- 
sion and time-span of the Planck data that was used for the real 
SZ cluster analysis of the previous section. 

The simulations were then analysed using the PwS algorithm 
in exactly the same manner as was applied during the analysis of 
the real data. For each of the 1 1 clusters, ten simulations were 
run. The dust template based on the Planck 857 GHz map was 
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the same for each of these ten simulations, but the CMB and 
noise realisations were different. 

The results of the simulations are shown in Fig. 7. In each 
panel, input parameters are indicated with a star and the re- 
covered parameter constraints from the ten different simula- 
tions are indicated by ten sets of differently coloured contours. 
Comparing with the input parameters, the recovered constraints 
for each of our clusters are clearly distributed about the input 
model and there is no indication of any significant bias due to 
dust contamination, noise bias, the Planck beams or the PwS ex- 
traction technique employed. 

6.2. AMI simulations 

As mentioned in Section 4, contamination from radio point 
sources is a significant issue at AMI frequencies (~16 GHz). 
Although the AMI LA observations are used to accurately find 
and model sources in the AMI SA observations, there is the pos- 
sibility of contamination from source residuals if this modelling 
is not perfect. 

In a similar manner to the Planck simulations described in 
the previous subsection, we have investigated potential issues 
associated with either residual foreground radio sources or with 
the AMI data-analysis methodology and instrument response us- 
ing simulations. The simulated input clusters were the same as 
used for the Planck simulations. 

To simulate the interferometric AMI observations, we used 
the in-house simulation package, profile (Grainge et al. 2002) to 
create the mock visibilities. In addition to the cluster signal, the 
simulations included primordial CMB fluctuations and Gaussian 
noise, the amplitude of which was chosen to match that mea- 
sured from the real observations. The simulation package also 
mimics the actual uv coverage and synthesised beam of the real 
observations. The point sources in each cluster were simulated 
using the best-fitting values from the analysis of the real data. 
These simulated observations were then analysed in the exact 
same way as was used for the real data. Once again, for each of 
the 1 1 clusters, ten simulations were performed. Here, the point 
source environment was kept the same for these ten simulations 
but the CMB and noise realisations were again different. 

The results of the simulations are shown in Fig. 8. In each 
panel, input parameters are indicated with a star and the recov- 
ered parameter constraints from the ten different simulations are 
indicated by ten sets of differently coloured contours. As was 
the case with the Planck simulations, the recovered constraints 
for each of our clusters are clearly distributed about the input 
model. Once again, there is no indication of any significant bias 
due to residual point source contamination, noise bias, the AMI 
uv coverage and resolution effects, or the extraction technique 
employed. 



7. Adopting individual pressure profiles as 
measured from X-ray observations 

The simulations presented in the previous section indicate that 
the discrepancies seen in the analysis of the real data cannot 
be easily explained by astrophysical contamination, instrumen- 
tal effects, or any issues associated with the SZ signal extraction 
techniques. It is then interesting to ask whether the discrepancies 
observed might be associated with the way in which the clusters 
have been modeled using the universal pressure profile (Amaud 
et al. 2010). 



For a number of clusters in the sample, we have high-quality 
determinations of the clusters' individual pressure profiles, as 
estimated from X-ray observations. Rather than adopting the 
Amaud et al. (2010) profile (which is essentially an average pro- 
file taken over many clusters), one might expect to achieve better 
consistency on a case-by-case basis if we use these individual 
best-fitting X-ray derived profiles in the SZ analysis. 

We have performed such a re-analysis for five clusters in our 
sample for which we have high-quality measured X-ray profiles. 
The clusters concerned are A1413, A1914, A2034, A2218, and 
A773. We fitted a GNFW pressure profile to the measured X- 
ray profiles and the results are presented in Table 7. We then 
re-analysed the Planck and AMI SZ data for these five clus- 
ters using the best fitting values of the profile shape parameters 
(c5oo,a,jS, y) as given in Table 7. The resulting constraints from 
this re-analysis are shown in Fig. 9. 

Comparing with the corresponding constraints for these five 
clusters in our original analysis, we see that the updated con- 
straints for A2034 have tightened significantly. This appears to 
be due to the fact that the previously used Amaud et al. (2010) 
GNFW profile was not a good match to this particular cluster's 
pressure profile, particularly in the central region of the clus- 
ter, where AMI is sensitive. The GNFW profile variant used to 
produce the updated constraints is a much better match to the 
measured X-ray profile and so the AMI data are better able to 
constrain the remaining cluster parameters. 

Apart from this single case, comparing with our original re- 
sults, there does not appear to be a systematic improvement in 
the agreement between the Planck and AMI constraints when 
we move from the Arnaud et al. (2010) profile to the best-fitting 
GNFW profile as measured from the individual X-ray obser- 
vations. This, and similar reasoning based on an adaptation of 
these modified profiles to the other clusters in our sample, sug- 
gests that a more significant widening of the parameter space 
describing the cluster profiles will be required in order to simul- 
taneously fit both the Planck and AMI SZ measurements for the 
entire cluster sample considered in this paper 



8. Conclusions 

We have studied the F500 - 6500 degeneracy from the SZ effect 
for a sample of 1 1 clusters (0.1 1 < z < 0.55) observed with both 
Planck and AMI. This is motivated by the fact that such a study 
can potentially break the well-known F-size degeneracy which 
commonly results from SZ experiments with limited resolution. 
Modelling the radial pressure distribution in each cluster using a 
universal GNFW profile, we have shown that there is significant 
overlap in the 2D posterior distributions for eight of the clusters. 
However, overall, AMI finds the SZ signal to be smaller in angu- 
lar extent and fainter than Planck finds. The derived parameter 
degeneracies are significantly different for the two instruments. 
Hence, where the constraints from the two instruments are mu- 
tually consistent, their combination can be powerful in terms of 
reducing the parameter uncertainties. Significant discrepancies 
are found between the Planck and AMI parameter constraints 
for the remaining three clusters in our sample. 

We have investigated the origin of these discrepancies by 
carrying out a detailed analysis of a series of simulations as- 
sessing the potential impact of diffuse thermal emission from 
dust and residual contamination from imperfectly modeled ra- 
dio point sources. Our simulations also include a number of sys- 
tematic effects associated with the two instruments in addition 
to primordial CMB fluctuations and thermal noise. We find that 
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Fig. 7. Recovery of SZ cluster parameters from simulated Planck observations (see Section 6.1 for details). Each set of recovered 
parameter constraints (contours with different colours) represents a different realisation of the instrument noise and primordial CMB 
fluctuations and the star shows the input parameter values. The inner and outer contours in each set indicate the areas enclosing 68 % 
and 95 % of the probability distribution. Any bias in the recovery of the input parameters averaged over realisations is negligible 
compared to the random errors. In cases where the contours do not close at the lower ends of the parameter ranges, the corresponding 
constraints represent upper limits only. 
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Fig. 8. Recovery of SZ cluster parameters from the simulated AMI observations in the presence of residual point source contam- 
ination from imperfectly modeled radio sources for each cluster in the sample. Also included in the simulations are the AMI uv 
coverage, the instrument beams and realisations of the instrument noise and primordial CMB. The different sets of contours indicate 
different CMB and noise realisations and the star shows the input parameters used to generate the simulated cluster The inner and 
outer contours in each set indicate the areas enclosing 68 % and 95 % of the probability distribution. 
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Table 7. Best-fitting GNFW shape and concentration parameters (cf. equation 3) derived by fitting the parameterised GNFW profile 
to the measured X-ray pressure profiles of five clusters in our 11 cluster sample. Note that p - 5.49 is fixed and a prior of y > is 
imposed. This latter constraint is enforced to avoid unphysical pressure gradients being allowed by the GNFW parameterisation. 
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Fig. 9. Constraints obtained from the re-analysis of the Planck and AMI observations for the five clusters for which high-quality 
X-ray observations are available. These re-analyses adopted the GNFW shape parameters which best fit the X-ray data as given in 
Table 7. Comparison with the con^esponding panels in Fig. 5 reveals no obvious improvement in the level of agreement between the 
Planck and AMI constraints. The X-ray sizes indicated here are also derived from the GNFW fits to the high-quality X-ray data. 
These are slightly different from the X-ray sizes plotted in Fig. 5, which were taken from the ESZ catalogue. 



the results of the simulations of both the Planck and AMI anal- 
yses are unbiased, confirming the accuracy of the two analysis 
pipelines and their corresponding methodologies. 

We have attempted to reconcile some of the discrepancies 
seen by re-analysing the Planck and AMI data adopting in- 
dividual best-fitting pressure profiles, as measured from high- 
quality X-ray observations for five of the clusters in the sample. 
However, we do not observe a systematic improvement in the 
agreement between the Planck and AMI parameter constraints 
when we perform this re-analysis. 

We conclude that: either (i) there remain unaccounted for 
systematic effects in one or both of the data sets beyond what are 



included in our simulations; or (ii) a further expansion of the pa- 
rameter space used to model the SZ cluster signal is required to 
simultaneously fit the Planck and AMI SZ data. Such further ex- 
pansion of the model parameter space, which we leave for future 
studies, could potentially include using the Planck and AMI data 
in conjunction with X-ray observations to find a global fit for the 
GNFW shape parameters, going beyond the GNFW parameter- 
isation to investigate other cluster profiles, and/or dropping the 
assumption of spherical symmetry for the SZ (and X-ray) emis- 
sion. 
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